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Abstract 

We study excitation and ionization in the 1.5 MeV/amu 0®'’'-Li collision system, which was 
the subject of a recent reaction-microscope-type experiment [Fischer et al, Phys. Rev. Lett. 
109 , 113202 (2012)]. Starting from an independent-electron model based on determinantal wave 
functions and using single-electron basis generator method and continuum distorted-wave with 
eikonal initial-state calculations we show that pure single ionization of a lithium iC-shell electron 
is too weak a process to explain the measured single differential cross section. Rather, our analysis 
suggests that two-electron excitation-ionization processes occur and have to be taken into account 
when comparing with the data. Good agreement is obtained only if we replace the independent- 
electron calculation by an independent-event model for one of the excitation-ionization processes 
and also take a shake-off process into account. 
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I. INTRODUCTION 


The investigation and differentiation of single and mnltiple electron processes has been 
a topic of considerable interest in atomic collision physics over many years. For one thing, 
studies in this area shed light on questions relevant for applied research, e.g., in radiation 
therapy, where the damage induced by swift ions is different in single- and multiple-ionization 
events ll). Multiple ionization of, e.g., water molecules results in fragmentation practically 
with certainty, while there is a high chance that the H 20 ’'' ion created after single-electron 
removal stays intact j2|. 

Much of the experimental and theoretical activity in the investigation of single and mul¬ 
tiple processes is, however, fueled by the fundamental interest in the few-body quantum 
dynamics at play. The seemingly simplest situation that can occur in an ion-atom collision 
corresponds to single-electron removal, in which one target electron is either captured by 
the (bare) projectile or promoted to a continuum state. However, single-electron removal is 
not necessarily a pure one-electron process: Another target electron may participate in the 
dynamics and end up in an excited state after the collision. In the case of capture, these 
so-called transfer plus target-excitation (TTE) events have been identified in cold target 
recoil ion momentum spectroscopy (COLTRIMS) and reaction microscope (ReMi) experi¬ 
ment^ which give access to the Q-value of a given reaction, i.e., the electronic energy loss or 
gain Excitation-ionization (El) was measured some time ago at very high projectile 

energy using (Auger-) electron spectroscopy Q. 

If the target atom is helium, which has often been the case in COLTRIMS and ReMi 
experiments, TTE and El are true two-electron processes. A somewhat different situation 
arises if more than one target shell is occupied before the collision. In this case, the single 
capture or single ionization of an inner-shell electron leaves the target ion behind in an 
excited state. In a naive independent-particle picture a two-electron process, in which an 
outer-shell electron is removed and an inner-shell electron is promoted to the state just 
vacated, has the same outcome. This type of two-electron TTE or El in principle becomes 
distinguishable from pure single inner-shell electron removal if the inner-shell electron is 


promoted to a target state that was vacant before the collision 3]- In practice, however, 
the resolution achievable in Q-value measurements may make it difficult to separate these 
two-electron processes from pure one-electron removal. Theoretical calculations are then 
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required for a full understanding of the situation. 

A recent joint experimental-theoretical work on 1.5 MeV/amu O®"*"collisions from lithium 
atoms was concerned with this problem 8|. The experiments were performed with the 
newly-developed ’MOTReMi’ apparatus, which combines a magneto-optical trap (MOT) 
to cool the lithium atoms with a ReMi to measure the reaction products. The recorded 
Q-value spectra for single ionization exhibit two distinct peaks, which can be associated 
with the ejection of the 2s valence electron and the ejection of an inner A'-shell electron, 
respectively. A continuum distorted-wave with eikonal initial-state (CDW-EIS) calculation 
for Li(2s) ionization was found to be in reasonable agreement with the measured electron- 
energy single differential cross section (SDCS) corresponding to the hrst peak, but a CDW- 
EIS calculation for Li(ls) ionization differed markedly from the data in the other channel. 
It was concluded that two-electron El, which was not taken into account in the calculation, 
is needed to explain the measurement. 

In this paper, we provide a theoretical analysis of the 0®^-Li collision system based on 
the independent-electron (lEL) model to scrutinize this interpretation. The lEL model is 
discussed in Sec. m In Sec. m we present lEL SDCSs for various pure ionization and 
El processes and compare them with the experimental data. Extensions of the lEL model, 
among them an independent-event (lEV) model for one El process, are considered to account 
for the quite substantial discrepancies. A summarizing discussion is provided in Sec. IIVI 
Atomic units characterized hy h = rUe = e = dvr^o = 1 are used unless otherwise stated. 


II. INDEPENDENT-ELECTRON TREATMENT OF PURE SINGLE IONIZATION 
AND EXCITATION-IONIZATION 

Within the semiclassical approximation (SCA) the 0®^-Li collision system is described 
by a time-dependent Schrodinger equation (TDSE) for the electronic Hamiltonian 

N N 

HM = j2if, + vr(t)]+J2%. ( 1 ) 

i=l i<j 

which consists of single-electron kinetic energy operators Tj, electron-nucleus interactions 
(which depend on time due to the SCA assumption of classically moving nuclei), and 
the electron-electron Coulomb interactions Wij. To the best of our knowledge an explicit 
solution of this correlated iV=3-electron problem has not been attempted yet. 
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The lEL model consists in replacing the Hamiltonian ([I]) by a one-body operator 


N 

He{t) —)■ hi{t) 

i=l 


( 2 ) 


snch that the TDSE separates into a set of single-particle eqnations for the three electrons. 
We assnme the single-particle Hamiltonian to be of the form 




( 3 ) 


where Zp is the charge nnmber of the bare projectile ion, and and Vp are the position 
vectors of the electron with respect to the target and the projectile center, respectively. They 
are related according to — R(t) with R(t) being the classical straight-line trajectory 

of the projectile relative to the target center. We note that the Laplace operator in Eq. ([3]) 
is taken with respect to the center-of-mass reference frame. The effective potential Vu 
represents the interactions in the (ls^2s) gronnd-state confignration of the lithinm atom. It is 
obtained from the exchange-only version of the optimized potential method (0PM) of density 
fnnctional theory j^, i.e., it inclndes electron-nnclens Conlomb interactions, screening, and 
exchange terms exactly and exhibits the correct asymptotic — l/r^ behavior, bnt it neglects 
electron correlations. 


We have solved the single-particle eqnations for the Hamiltonian (|3]) and the initially 
occnpied Li(ls) and Li(2s) orbitals nsing the two-center basis generator method (TC-BGM) 


with a basis that consists of the Is — 4/ target states, the Is — 4/ 
states, as well as 71 BGM psendo states to acconnt for ionization 


hydrogenlike) projectile 


10| . All basis states are 


endowed with electron translation factors to ensnre Galilean invariance. Resnlts for the total 
excitation (Pi^( 2 s)), captnre (Pi 5 ^ 2 s)) ionization (pj°J( 2 s)) probabilities are shown in Fig. [1] 
as fnnctions of the impact parameter b. The probabilities are calcnlated as follows: 
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Pls^fk = lifkl'lplsitf))]^ ( 4 ) 

P2s^h = l(/fc|^/’2«(t/))r (5) 

p17 = Y1 

fkSTJkj^ls 

P’S = E P2„* (7) 

/feeT,/fe 7 ^{ls ,24 

pT 7 = Pis^h ( 8 ) 

fk€P 

pT7 = Y 

fk&p 

Pi 7 = ^-Y~ Y 

fk^'^ fk^P 

P 2 T = ^-Y “ Y (11) 

fk^P fk^P 


where |' 0 is( 2 s)(V)) denote the solutions of the single-particle equations corresponding to the 
Is(2s) initial states and propagated to a sufficiently large hnal time tf after the collision, 
and \fk) € T(P) are hnal target (projectile) states. Note that in Eqs. (fTOj) and flTT]) we use 
the unitarity of the TC-BGM solutions of the single-particle equations. 

As expected at a projectile energy as high as 1.5 MeV/amu, capture is very weak except 
for the Is initial state in close collisions. It will be neglected in the following, i.e., we will 
identify electron removal from the lithium atom with ionization into the continuum. We 
further observe in Fig. [1] that ionization strongly dominates for the case of the Is initial 
state, while for the Li(2s) initial state excitation takes over in distant collisions. The details 
of our calculations show that excitation to 2p is the strongest channel. We also include in 
jjg. [1] ionization probabilities obtained from a CDW-EIS calculation for the Hamiltonian ([3]) 
111 . Il2| . They are in very good agreement with the TC-BGM results for the case of the Is 
initial state, and still in acceptable agreement for 2s—keeping in mind that the CDW-EIS 
method is perturbative in nature and the perturbation parameter r] = Zp/v is close to one 
{rj = 1.03 for the projectile speed v = 7.75 a.u.). 


In order to make contact with the experiment [8] and take El processes into account, 
we have to reinstate many-body aspects of the collision system. Consistent with the lEL 
model this is done by assembling the solutions of the single-particle equations in the form of a 
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Slater determinant. If the three-electron hnal states are also taken to be Slater determinants, 
the transition probabilities of interest can be obtained without further approximation from 
combinations of determinants constructed from one-particle density matrix elements 1^ 


N 




( 12 ) 


i=l 


In this work, we are interested in processes, in which exactly one vacancy is created in the 
lithium atom. They correspond to probabilities for hnding two bound single-particle target 
states, say |/i) and I/ 2 ), occupied and all the others vacant, and can be calculated according 
to |l3| 


P 


Efc fk — 


/1/2 


= P. 


/1/2 “ X] Phf2fk- 

fk&T 


(13) 


In the expression on the right hand side of Eq. (fT3l) P/ 1/2 denotes the inclusive probability 
for hnding two electrons in the subconhguration \ f_ 1 f 2 ) while nothing is known about the 


hnal state of the third electron. As shown in Ref. 


13j this inclusive probability is given as 


the determinant constructed from the 2x2 density matrix corresponding to 1 / 1 / 2 )- Pfi/ 2 /fc is 
the exclusive probability to hnd the three electrons in the completely specihed conhguration 
1 / 1 / 2 /fc) and is given as the determinant of the 3x3 density matrix corresponding to |/i/ 2 /fc)- 
Since the sum in Eq. ffT3D runs over all bound target states the diherence of both terms 
corresponds to the statement that one of the electrons has been removed from the lithium 
atom. 

For the explicit evaluation of Eq. flTSD we have to take the spin projections (/, /) of the 
electrons into account and consider spin orbitals that we denote, e.g., by writing \fk /). 
Since the Hamltonian (j3]) is spin independent we have {fk t l/’i t) = {fk i iV’i /) = (/fcl'^i) 
and {fk t I'f’i /) = {fk i I'fi t) = 0. We choose the initial state of the lithium atom to be a 
Slater determinant built from the (Is/ Is/ 2s/) spin orbitals. 

With these preparations we are ready to consider the probabilities that correspond to the 
experimentally distinguishable processes 8|]: 

(i) 2s vacancy production in which the Li’*' ion is found in its (Is^) ground state: 


pvac _ pXIfe At 


(14) 


(ii) Is vacancy production Pff^, in which the Li’*" ion is left in the excited conhguration 
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(Isnl) with n > 2: 


pvac _ \ ^ /pZlfe/fc'l-i pZlfe /fet I pZlfe/fet\ M c;'! 

^is - 2^ )■ 

/iSTJi^ls 

All terms on the right hand sides of Eqs. flT^ and flT^ can be computed using the 
prescription flT^ . It is instructive to work out the determinants analytically. With the 
dehnitions (jl]) to (fTTj) and 

pfr = Pls^ls (16) 

ptr = P2s^2s (17) 

one obtains 


p,T = Ar'+AT+( 18 ) 

Pi7" = Ar' + Peii + Pei 2 + Pi7 + APi^*' (19) 


with 


pexcl _ / elast'i2 ion 

-'2s ~ Inis ) P2s 

pex _ „elast ion 
^2s ~ Pis Pis n2s-)-ls 

pexcl _ 9 elast ion elast 

-'^Is ~ ^Pls Pis P2s 

t^EIl - 2p^^ P2s 

p _ o^elast^exc^ion 

t^EI2 — ^Pls Pis P2s 

Pl 7 = 2 pSV“P 2 s^ 1 s 

and correction terms and stemming from off-diagonal elements of the density 

matrix that reflect the antisymmetry of the three-electron wave functions. These terms as 
well as a few details regarding the derivation of Eqs. (fT8|) and (IT^ are given in the Appendix. 

Equations fl20|) to fl2^ have straightforward interpretations and can also be obtained from 
a simple multinomial analysis of the three-electron problem that is based on associating each 
electron with a single-particle probability. Expressions (12(11) and (|2^ correspond to exclusive 
ionization events, in which the non-ionized electrons remain bound in their initial states. The 
probabilities fl23l) and fl2Tl) correspond to two-electron El processes, in which one P-shell 
electron does not change its state, while the other two electrons are either excited or ionized. 
Finally, Eqs. (|2T1) and fl2^ describe exchange processes which arise as a consequence of 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


7 


the indistinguishability of the electrons. In practice, these exchange terms as well as the 


correction terms and are negligible. 


and APi'^^^ 

This is demonstrated in Figs. [2] and |3l Fignre |2] displays the 6-weighted probability 


for 2s vacancy prodnction according to Eq. o. Differences to exclnsive 2s ionization 
according to Eq. (120|) are too small to be visible mainly becanse the exchange amplitnde for 
a transition from the 2s to the Is state in the correction terms is very small. We have also 
inclnded in Fig. [2l This probability is interpreted as the single-ionization probability 
in a one-active-electron (OAE) model in which the two iF-shell electrons are assnmed to 
be frozen thronghont the collision. Obvionsly, the OAE resnlt differs from the other cnrves 
only at relatively small b, where the Is electrons do nndergo transitions with non-negligible 
probabilities. 

In Fig. [3] we show 6-weighted probabilities for processes that resnlt in the creation of 
one iF-shell vacancy. Peii [Eq- (ES])] is almost as strong as Pi^^^ [Eq. (122]) ]. whereas Pei 2 
[Eq. fl2^ ]. which involves the excitation of a iF-shell electron, is considerably weaker. The 
snm of these three probabilities is in almost perfect agreement with the fnll resnlt for Is 
vacancy prodnction according to Eq. (na — demonstrating that also for Is ionization the 
exchange and antisymmetry correction terms in Eq. flT^ can be neglected. 

In many calcnlations for helinm targets the qnantity 2p\°“ was nsed to caicniate ionization 
cross sections (see, e.g.. Ref. j^, 15[). This procednre again corresponds to an OAE model, 
in which the factor of two arises since one does not know which of the two iF-shell electrons 
is active and which is passive. We have inclnded a 2p^°^ cnrve in Fig. [3l It is seen to be 
mnch larger than the other probabilities, which indicates that the assnmption of jnst one 
active iF-shell electron is not realistic for the 1.5 MeV 0®+-Li collision system. The main 
reason for this is that it is very nnlikely for the initial 2s electron not to nndergo a transition 
in a relatively close collision. Pnt another way, if the assnmption of two passive electrons is 
nnrealistic, the qnantity 2^^°^^ is contaminated by mnltielectron processes. 


III. DIFFERENTIAL CROSS SECTION RESULTS AND EXTENSIONS OF THE 
lEL MODEL 


Having analyzed the relative strengths of varions ionization processes within the lEL 


model we now tnrn to the SDCSs for the two reaction channels discnssed in Ref. j^. Ideally, 






we would calculate them on the basis of Eqs. flT^ and (IT^ and the TC-BGM single-particle 
solutions. However, this is not possible, since it is difficult to extract electron-emission- 
energy (Eg) differential information from the population of the BGM pseudo states. Instead, 
we use the GDW-EIS method to calculate single-particle ionization probabilities Eg) 

for i = Is, 2s and combine them with the 6-dependent TG-BGM probabilities for excitations 
and elastic transitions according to the simplihed Eqs. (I20l) . (|22|) . fl23|) . and fl2^ . The 
resulting 6 and Eg dependent probabilities are integrated over the impact parameter in 
order to obtain the Eg-differential SDGSs. Given that (i) the exchange and antisymmetry 
correction terms in Eqs. fllSp and (IT^ were shown to be small, and (ii) the GDW-EIS total 
ionization probabilities ^““^(6) were shown to be close to their TG-BGM counterparts this 
procedure is unlikely to produce errors in addition to those stemming from the limitations 
of the lEL model. 

Results are shown in Fig. H] in comparison with the experimental data of Ref. jsl. For 
the 2s channel the OAE and exclusive ionization SDGSs are very similar (cf. Fig. [2]) and in 
reasonable agreement with the data, which were deemed low at energies Eg > 20 eV due to 
acceptance limitations of the spectrometer 8j. We note that the present results are in very 


f-fl. 


good agreement with the GDW-EIS calculations shown in Fig. 2(b) of Ref. |8i, which we did 
not include in our hgure for the sake of clarity. 

For the Is channel the situation is more involved. Our OAE result agrees well with the 
GDW-EIS calculation reported in Ref. 8j (not shown), but it differs in both magnitude 
and shape from the experimental data. As mentioned earlier, the OAE cross section is con¬ 
taminated by multielectron processes, particularly by multiple ionization, which is excluded 
in the measured coincidences of electrons and singly-charged ions. Hence, this quantity is 
not very useful for the analysis of the experimental data. On the other hand, exclusive Is 
ionization is much smaller than the measurements at all electron energies. The same is true 
for the El process described by Eq. (I23l) . which consists of Is ionization and 2s excitation. 
Similarly to exclusive Is ionization this El process reflects the energy dependence of the 
OAE Is ionization curve, which can be inferred from Eqs. (1221) and (|23|) . 

Likewise, the other El process [Eq. (I2T)) ] which involves 2s ionization and Is excitation 
reflects the energy dependence of the 2s ionization curve. If one adds it to exclusive Is 
ionization and Is ionization with 2s excitation to calculate the total lEL SDGS for Is- 
vacancy production one obtains a curve whose shape has some similarity with the energy 
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dependence of the experimental data, but lies significantly below them except at the highest 
electron energies. We thus have to conclude that the lEL model is not sufficient to explain 
the measurements. 

In an attempt to understand the discrepancy without embarking on a full solution of 
the correlated three-electron problem we have considered the following scenarios. First, one 
may argue that for the El process corresponding to Eq. (|2^ a sequential lEV picture of 
the collision dynamics is more realistic than the lEL model. In such a scenario the outer 
(2s) electron is removed first with the result that the inner (Is) electrons rearrange before 
one of them is promoted to an excited state of the Li’*' ion. Accordingly, we also considered 
a Hamiltonian in which the 0PM potential Wi is replaced by the 0PM potential Vu+ for 
the (Is^) conhguration which decays like —‘Ijri asymptotically. The Is-excitation with 2s- 
ionization SDCS in the lEV model was then calculated according to 

da 


dEp 


= 2 / d^bpf:f(b)pzi(b)p‘i‘(b,E,), 


( 26 ) 


lEV 


where pf^^(b) and pfg^(b) denote single-particle elastic and excitation probabilities obtained 
from TC-BGM collision calculations for 0®’''-Li’'', and p^ 2 s^(b, Ep.) is the electron-energy differ¬ 
ential CDW-EIS ionization probability for 0®’''-Li that we also used in the lEL calculations. 
We note that we did not consider the lEV counterpart of the El process fl2^ . in which a 
A'-shell electron is removed and the 2s electron excited, since it does not seem plausible 
to assume that the wave function of the outer 2s electron adjusts to an ionic (ls2s) con¬ 
figuration after inner-shell ionization. Besides, we know from the analysis presented above 
that such an lEV probability would result in the same Eg-dependence as OAE Is ionization, 
which is at odds with the measured data. 

Secondly, our TC-BGM basis only includes target states up to principal quantum number 
n = A. It is impossible to push the calculation further without severely restricting the repre¬ 
sentation of the continuum. In order to estimate the contribution from higher excitations we 
have assumed that they scale like 1/n^ which allows for extrapolating them to n —)■ oo 16|. 

Finally, we have estimated the contribution from a correlated two-electron process, in 
which a Is electron is excited and the 2s electron is emitted due to shake-off. We have calcu¬ 
lated the (Eg-differential) shake probability by projecting the Li(2s) state onto continuum 
states of the single-center Hamiltonian 


k 


shake 


= -iA 


A + i[l - (1 + Z,.r,)e-2Z„.| 

n n 


(27) 
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For Zt = S the total potential of this Hamiltonian corresponds to the sum of the Coulomb 
potential of the lithium nucleus and an electrostatic potential due to the presence of one 
hydrogenlike Is electron, i.e., the effects due to the presence of the excited electron and the 
repulsion between the Is and 2s electrons are neglected. We have repeated the shake-off 
calculation with continuum states of the Li’*' Hamiltonian used in the lEV calculation and 
have found very similar results indicating that these simplifications are of minor importance. 

Figure [5] displays the results of these extensions. The lEV SDCS (j26l) gives the strongest 
contribution. The correction due to the n ^ oo extrapolation amounts to about 20% and 
is included in the cross section curve shown. The strongest partial channels are those corre¬ 
sponding to the dipole-allowed excitations into 2p (~ 44%) and 3p (~ 13%). Interestingly, 
the lEV SDCS is much stronger than its lEL counterpart for the process (|2T|) shown in 
Fig. m This mirrors the fact that the single-particle excitation probabilities for the lithium 
ion (pis+) are larger than those for the atom (pi“), which can be understood by inspecting 
the orbital densities of the single-particle states involved. As shown in Fig. Elthe Li''"(ls) and 
Li(ls) densities are practically indistinguishable on a linear plot. They only differ in their 
asymptotic decays which are determined by the energy eigenvalues (indicated as horizontal 
lines in Fig. [6]). By contrast, the excited states are markedly different with those of Li'*' 
being more compact than those of Li and thus easier accessible for an initial iF-shell elec¬ 
tron. Figure [6] shows that in the case of the neutral lithium atom the Is and 2p states have 
only very little spatial overlap. These observations favor the lEV picture, because the ionic 
excited states are more realistic final states for El than the neutral atom states. Figure [6] 
also displays the Li and Li'*' 0PM potentials. The shoulder in the neutral atom potential is 
a reflection of the shell structure. 

Coming back to Fig. |5l we observe that the Is ionization with 2s excitation SDCS cor¬ 
responding to Eq. (|23|) increases by a little less than 20% when the contribution from exci¬ 
tations to n > 5 is included, and becomes very similar to the exclusive Is ionization curve 
(which is repeated in Fig. [5]). Again the 2p excitation channel is the strongest one, con¬ 
tributing almost half of the cross section. The shake-off SDCS, obtained by replacing 
in Eq. (|2T)) by the shake probability, applying the n ^ oo correction to and integrating 
over the impact parameter, falls off steeply with increasing electron energy, but it does give 
a non-negligible contribution at low E^. When added to the three other contributions shown 
in Fig. [5]we obtain the thick full curve. The agreement with the measurements is not perfect. 
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but it is very good down to Ee'^2 eV. 


IV. SUMMARY 


In this work, we have used the independent-electron model to analyze the contributions 
from single ionization with and without additional target excitation to the recent MOTReMi 
SDCS measurements for 1.5 MeV/amu 0®’''-Li collisions. Our results, obtained by combining 
TC-BGM and CDW-EIS single-particle probabilities, largely confirm what was speculated in 
Ref. jsj: Two-electron excitation-ionization processes occur and must be taken into account 
when comparing with the measurements for the Is channel. However, the independent- 
electron model is not sufficient to explain the data. Rather, we had to replace it by an 
independent-event calculation for one of the El processes, i.e., we had to take relaxation 
effects into account, and we also had to include a shake-off process to obtain good agreement 
with the measurements. Accordingly, we have to conclude that electron correlation effects 
play some role in this collision system. 

It would be of interest to extend the experimental and theoretical studies by varying the 
collision parameters projectile charge and velocity. This would provide additional insight into 
the limitations of the independent-electron model and the role of relaxation and shake-off 
processes. Ultimately, the theoretical challenge consists in carrying out a full three-electron 
calculation for this collision system. 
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Appendix 

According to Eqs. flTOl) and flTTl) the 2s vacancy production is given as 

p^T = n.t.4 - E n.n4/.t- (A.i) 

h&T 
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The inclusive probability is the determinant of the 2x2 density matrix corresponding 

to the conhguration (Is t Is i): 


_ (Is t |7H^/)|1 s t) (1st |7^(t/)|ls;) 

-tlstls4, - ) 

(Is t |7H^/)|1 s t) (Is t 1 ) 

where [cf. Eq. flT^ ] 

fitf) = \^lst{tf)){^lstitf)\ + + |V’2«t(^/))(V’2st(^/)|- 


(A.2) 


(A.3) 


Inserting flA.3j) into flA.2p . exploiting spin orthogonality and omitting the time argument for 
convenience we obtain 


P^ 


Isflst 


\{ls\^Pu)\^ + \{ls\^l^,s)\^ 0 

0 l(ls|V'i.)P 


= {p 


,elast N 2 I ^elast ^ 

Is ) +Pls P2s^ls- 


(A.4) 


In the last step we have used the dehnitions dlj), ([S]), and flTBD . 

The exclusive probability Pisfuiftt determinant of the 3x3 density matrix corre¬ 

sponding to the conhguration (Is t Is t /fc t)- Using similar arguments as above we can 
write 


\{ls\'ljjis)\'^ + \{ls\^p 2 s)\‘^ 0 (ls|V’ls)(V’l^|/fc) + (ls|^/>2s)(^2s|/fc) 

UlstWfct = 0 |(ls|'01s)P 0 

(/fc|^/’ls)(V’ls|ls) + (/fc|^ 2 s)(V’ 2 s|ls) 0 \{fk\'lpls)\‘^ + |(/fc|^ 2 s)P 

= {pfr^YP2s^f^ + pfr^[pis^f,,P2s^ls - ((ls|'01s)(V’ls|/fc)(/fc|^/’2s)(V’2s|ls) + C.C.)], (A.5) 


where c.c. denotes the complex conjugate. One easily infers from the last expression of 
Eq. (IA.5D that ~ ^ l/fc) ~ 1^®) required by the Pauli principle. 

Inserting flA.4p and flA.51) in (lA.ip we obtain 



+ pfs^^ 5^((ls|?/'ls)(V’ls|/fc)(/fc|^2s)(V’2s|ls) +C.C.). (A.6) 

fk(^T 

In the last step we use the dehnitions for single-particle ionization fllOp and fllip and remem¬ 
ber that capture is negligible to obtain Eq. (ITSP with the abbreviations (|20P and fl2Tp and 
the correction term 


APC*' = P?r‘ ({lsW'i.)«'i,IA)(Alfc)(fe|ls) + C.C.). (A.7) 

fk&T 


13 




















Equation ffT^ for the Is vacancy production can be derived using similar arguments. 
The calculation is straightforward but lengthy. One Ends for the correction term: 

^panti ^ _pion ^ ((Is(V'l. |/fc) (/fc | V’2s) (V'2. | Is) + C.C.) 

fkeTJkj^ls 

+ PlT '^i{^s\^pu){'lpls\fk){fk\'4^2s){i^2s\^s) + C.C.) 
fk<^T 

+prE E {{fk\i^ls){i^ls\fl){fl\i^2s){i^2s\fk) +C.C.). (A.8) 
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FIG. 1. (Color online) Total excitation, capture, and ionization probabilities according to Eqs. ([6]) 
to (fTTI) for 1.5 MeV/amu 0®^-Li collisions obtained from TC-BGM calculations for Li(ls) (upper 
panel) and Li(2s) (lower panel) initial states. CDW-EIS ionization probablilities are displayed in 
addition as (black) dotted curves. 


16 













FIG. 2. (Color online) Impact-parameter-weighted probabilities for 2s vacancy production 
[Eq. (fni) ]. exclusive 2s ionization [Eq. (I20p ]. and single-particle 2s ionization [Eq. dlip ] for 1.5 
MeV/amu 0®^-Li collisions. All curves shown are obtained from TC-BGM calculations. The 
curves displaying the probabilities (I14p and (12011 sit on top of each other. 
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FIG. 3. (Color online) Impact-parameter-weighted probabilities for Is vacancy production 
[Eq. exclusive Is ionization [Eq. (1221) ]. Is ionization with 2s excitation [Eq. (I23p ]. Is ex¬ 

citation with 2s ionization [Eq. ([24)) ]. the sum of (I22p . (l23l) and (|24l) (X))) single-particle Is 
ionization multiplied by two [Eq. (1101) x 2] for 1.5 MeV/amu 0®'’'-Li collisions. All curves shown are 
obtained from TC-BGM calculations. 
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FIG. 4. (Color online) Single differential cross sections (SDCSs) for single-particle 2s ioniza¬ 
tion [Eq. dill) ], exclusive 2s ionization [Eq. (I20p ]. single-particle Is ionization multiplied by two 
[Eq. (I10p x2]. exclusive Is ionization [Eq. (1221) ]. Is ionization with 2s excitation [Eq. (1231) ]. Is ex¬ 
citation with 2s ionization [Eq. (I24p ]. and the sum of (|22ll . (12311 and ([24) (^) as functions of the 
electron energy for 1.5 MeV/amu 0®"’“-Li collisions. Experimental data: 8[. 
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FIG. 5. (Color online) Single differential cross sections (SDCSs) for 2s ionization with ls+ excitation 
extrapolated to n —>■ oo [Eq. (I26p for all n], the shake-off model as explained in the text, exclusive 
Is ionization [Eq. (l22l) ]. and Is ionization with 2s excitation extrapolated to n —oo [Eq. ([23)1 for 
all n[ as functions of the electron energy for 1.5 MeV/amu 0®'’“-Li collisions. The full thick curve 
dubbed ’final total EE is obtained from adding the four other calculated SDCSs shown here and 
represents our final result for total EL Experimental data: 8|. 
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FIG. 6. (Color online) 0PM potentials and Is and 2p radial orbital densities for Li and Li^. The 
Li(ls) and Li+(ls) densities sit on top of each other. The horizontal lines indicate the binding 
energies of the Is and 2p orbitals. 
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